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ABSTRACT 

We give a new proof of the renormalizability of a class of matter field theo- 
ries on a space-time lattice; in particular we consider c/> 4 and massive Yukawa 
theories with Wilson fermions. We use the Polchinski approach to renor- 
malization, which is based on the Wilson flow equation; this approach is 
substantially simpler than the BPHZ method, applied to the lattice by Reisz. 

We discuss matter theories with staggered fermions. In particular we analyse 
a simple kind of staggered fermions with minimal doubling, using which we 
prove the renormalizability of a chiral sigma model with exact chiral symmetry 
on the lattice. 
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Perturbative quantum field theories on the lattice |I[ present a few difficulties that are 
absent in other regularization schemes; their renormalization has been studied by Reisz 
H, [|, £| in the BPHZ approach || || . Beyond the usual technical difficulties of the BPHZ 
approach, requiring an analysis of the topological structure of the Feynman graphs, there 
are those specific to lattice theories, due to the periodicity of the Green functions in 
momentum space. In the first place, it is non-trivial to prove a power-counting theorem 
on the lattice |2[ , since the Feynman integrands are not rational functions of the momenta, 
but of trigonometric functions of the momenta. In the second place one has to show that 
the subtraction procedure can be implemented with counterterms which are local on the 
lattice. 

In lattice gauge theories there is an additional difficulty, the presence of an infinite 
number of irrelevant interaction terms in the bare action; Reisz [J| has shown that these 
vertices do not modify the renormalized Green functions computed in the continuum. 

In theories containing fermions there is the doubling problem for fermions on the 
lattice J?], [|; introducing the Wilson term || in the fermionic action the doublers decouple 
in the continuum limit, but chiral symmetry is broken in a hard way by the Wilson term. 
Reisz's results apply to Wilson fermions p], [3|, 3, [TCfl ; these results have not been extended 
to fermionic models with doublers. While chiral invariance is not maintained by Wilson 
fermions, it is possible to construct models with doublers having chiral invariance on the 



lattice, for instance with the Kogut-Susskind |L1J] staggered fermions, or with simpler 
versions of staggered fermions [Tj|; it has not yet been shown that these models are 
renormalizable. 

A simple approach to renormalization in continuum quantum field theory has been 
initiated by Polchinski [Tl|, in the spirit of the Wilson renormalization group ]T3|; this 
approach has been further simplified and improved in JIB , |I"7j , [18| . One studies an effective 
action which, besides the ultraviolet cut-off A , has an 'infrared' cut-off A. The effective 
action satisfies perturbatively a linear differential equation in A, called the flow equation. 
One can impose mixed boundary conditions on the effective action, fixing the relevant 
terms of the effective action at a renormalization scale A = A 1; while the irrelevant terms 
are fixed at a scale A = Ao- 

In this paper we study the renormalizability of scalar and Yukawa theories on the 
lattice; we adapt the Polchinski approach to the lattice, introducing an 'infrared' cut-off 
A, which restricts the propagators on a periodic band in momentum space, at a distance 
of the order of A from the poles of the propagators. 

We extend Reisz's proof of renormalizability to a class of Yukawa theories with a 
simple kind of staggered fermions |~2], |i~3f . 

In the first section we study the massive 4 field theory on a hypercubic lattice; we 
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prove that it is multiplicatively renormalizable, using the flow equation for the amputated 
connected Green functions; we show that the addition of irrelevant terms to the action 
does not change the renormalized Green functions. 

In the second section we prove the renormalizability of a general class of matter field 
theories on a lattice including scalars, fermions and auxiliary fields. 

In the third section we apply these results to Wilson fermions and a simple kind of 
staggered fermions with minimal doubling; using the latter we prove the renormalizability 
of a chiral sigma model with explicit chiral invariance on a hypercubic lattice. The proof 
is obtained by rewriting the model on a reduced lattice, on which the results of the 
previous section apply. 



1 Renormalization of the massive </> 4 model on a hy- 
percubic lattice. 



1.1 Introduction. 

Consider an infinite space-time hypercubic lattice with lattice spacing a and sites in 
x ii — n fi a , where n M are integers, fJ> — l, ...,4. The bare action is 

„2 ™(0)2 _(0) 

s = \ E 4 0) (20i o) - <t>% - <t>%) + ^ E(V^° )2 + V^° )4) + s - (L1) 

where <f>^ is the bare field; rrS^ and are the bare parameters. S irr is an irrelevant 
term, which will be specified below; it will be chosen in such a way that the action 
preserves the hypercubic space-time symmetry and the discrete symmetry — > — 
The action reads, in momentum space, 

1 r n(°) r 4 4 

S =t ^) { - k ){k 2 + m^)^\k) + 9 — (27r) 4 4(^A; l )n0 {O) (^) + ^(1.2) 

/ Jk 4! Jfci,...,fc 4 1 i=1 



where k ^ = ^sin-j^, 5%{k) = j§^J2x e lkx is the periodic delta function, 

z d 4 k 
-f (2^) 4 



r — r a d 4 k 
Jk — J-S ( 2?r )4- 



The Fourier transform of <p x is <fi(k) = ^4 J2x z~ lkx 4>x-, and <p x = J k 4>(k)e lkx . 

The Fourier transform (f>(k) and the kinetic operator k are periodic under k^ — > 
fc M + ^ The momenta can be restricted to the Brillouin zone \k^\ < ^, provided one 
requires momentum conservation at the vertices only modulo — . 
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Ao = - (1.3) 
a 

is a kind of ultraviolet cut-off on the cartesian components of the momenta. 
In terms of the renormalized field <p x = Z~^(f)^ the action becomes 

S[4] = \j k 4>{-k){k 2 + m 2 )0(A;) + SM] (1.4) 

with 

-i 4 4 

= o / 0(-A:)(ci + c 2 P)0(A;) + ^ / (2vr) 4 4(^ A;,) <t>{h) + (1.5) 

z Jk 4! JkiMMM i i= i 

where 

ci = Zm (0 » 2 - m 2 ; c 2 = Z - 1 ; c 3 = # (0) Z 2 (1.6) 

In terms of the loop expansion parameter H the renormalization constants c» have the 
series expansion q = g5 ij3 + YlaLi'h 1 ■ The cf^'s are constants depending on A , m, g 
and on the renormalization conditions. A standard set of renormalization conditions is, 
in terms of the proper vertices, 

1 4 <9 2 r 

r 2 | P =o = ™ 2 ; | P =o = i ; r 4 | p=0 = (? (1.7) 

The parameters m and (7 are called renormalized parameters. A more general set of 
renormalization conditions will be specified later. 

S^,,^] contains irrelevant terms which are defined on the lattice; it will be chosen 
linear in a parameter y e [0, 1], Si rr [0] = ySi„\4>\\y=i- The dependence of Si rr [<f>} and of 
the Green functions on y will not be explicitly indicated in the following. A typical choice 
(though not the most general) is S irr [(f)\ = y J2x d2 n o 2n( Pl n , where d 2n , with n > 3, 
are independent parameters. A more general form for Si rr will be specified later, when 
we will choose the renormalization conditions. Obviously a trivial choice is Si„ = 0. 

Lattice perturbation theory is done using the Feynman rules, in which the propagator 

is 

Dip) = ^]— 2 (1.8) 

and the vertices are contained in Si. The dependence of D(p) on the ultraviolet cut-off 
A (1.3) is not explicitly indicated. 

The integrands of the Feynman integrals are periodic functions of the momenta; it is 
not possible to define the degree of divergence of a graph by considering its behavior for 
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very large momenta, so that the usual power-counting theorems on the continuum do not 
apply on the lattice. A detailed analysis of the structure of the lattice diagrams shows 
that it is nevertheless possible to define a lattice degree of divergence and to establish a 
power-counting theorem on the lattice 0. 

The fact that the <p A model is perturbatively renormalizable (multiplicatively) on 
the lattice means that, for n > 2, the n-point connected Green functions computed 
with these Feynman rules admit finite limit, order by order in h, in the continuum 
limit a — > 0, provided the constants c\ , i = 1,2,3 are chosen in such a way that the 
renormalization conditions are satisfied, and provided the irrelevant vertices contained 
in Si rr are chosen to satisfy a bound to be specified later. A proof of multiplicative 
renormalization of <p A using BPHZ has been given by Reisz @, where it was also 
shown that certain irrelevant terms, like Si rr [(j)} = J2 X d2n^ 2n <P x n , do n °t change the 
renormalized Green functions. Reisz has also shown that there is a considerable freedom 
in the choice of the lattice propagator 0, |j. We will return to this point at the end of 
this section. 

We will give another proof of these facts using an approach to renormalization started 
by Polchinski [|14| and improved in |16| . 

The idea followed in Jl6| to study renormalization consists in introducing a class of 
theories, labelled by a continuous parameter A 6 [0, Ao], differing from the previous 
theory only for the propagator, called D\. For A = the two propagators are equal, 
Dq(p) = D(p), while for A = A one has D^ = 0, so that for A = A perturbation theory 
is trivial, since the proper vertices coincide with the bare vertices. The Green functions 



of this class of theories satisfies a 'heat' flow equation, first considered by Wilson |L5| and 
further studied by Polchinski |TJ . 



Renormalization of the theory is studied investigating the 'heat' flow of the Green 
functions from A = A to A = 0. The theory at A = is the one in which one is really 
interested, that is the one with the action (1.4 — 5). 

Mixed boundary conditions are applied to determine this flow. The renormalization 
conditions on the relevant vertices (see e.g. the renormalization conditions (1.7)) are 
interpreted as boundary conditions at A = 0; the irrelevant vertices are fixed by the 
boundary conditions at A = A . 

Using the flow equations one can prove by induction inequalities on the Green func- 
tions. From these inequalities one can show that the Green functions are bounded, and 
finally that they converge, for Ao — > oo. 

In []TJ| one uses a smooth cut-off on the momenta, which vanishes below an infrared 
scale A and above an ultraviolet scale A . 

On the lattice the situation is somewhat different, since the Feynman integrands are 
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periodic functions of the momentum, with periodicity 2A (see (1.3) ). Some care must 
be taken in defining what is a 'small momentum' or a 'large momentum'. The propagator 
of a field is a periodic function on the Brillouin zone; roughly speaking, a momentum is 
small if it is close to a maximum of the propagator; it is large if it is close to a minimum 
of the propagator. The maximum of the propagator (1.8) for the scalar field is at p^ = 
in the Brillouin zone; due to periodicity it has infinite other maxima on R 4 , for instance 
in p M = 2A . On the other end the point p^ = A is a minimum of the scalar propagator. 
Then p^ = A is a 'large momentum', while p M = 2A is a 'small momentum'. The notion 
of region of large momenta is substituted on the lattice by the notion of periodic band 
at a distance of the order of A from the maxima of the propagators. 

In the case of the scalar propagator, such a band can be identified with a region 
around the boundary of the Brillouin zone. In the case of staggered fermions, which 
we will discuss in a later section, the propagator has a maximum in p M = and other 
maxima along the border of the Brillouin zone, so the band of large momenta does not 
coincide with the boundary of the Brillouin zone. 

We will show that there is a natural way to extend the flow equation on the lattice; 
the infrared cut-off A regulates the distance of the periodic band mentioned above from 
the maxima of the propagators. For A = the band coincides with the full Brillouin 
zone; as A increases, the band recedes from the maxima of the propagators of the fields, 
until for A = Ao it becomes vanishingly small. 

1.2 The flow equation. 

Introduce an 'infrared' cut-off A in the propagator, with < A < A , 

B,W = ^ (1-9) 



IS 



where, for A > 0, K\(p) is a C°° function of fa which vanishes for fa < -|, i 
monotonously increasing for ^| < fa < |f and is equal to 1 for fa > ||. For A = we 
define K = 1. 

K\(p) is a periodic function on the Brillouin zone; for A small, it cuts off a slightly 
squeezed ball around p = 0. For A close to A , it is non-vanishing only in a slightly 
squeezed ball around p = (-,...,-). 

This propagator satisfies the following conditions 

D Ao = ; £> (p) = _L^. (i.io) 

p l + m l 

D\ and K A depend on the ultraviolet cut-off A ; here and in the following the dependence 
on A will not be explicitly indicated. 
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Observe that lim Ao ^ooD A (p) is rotation ( 0(4) ) invariant. 
The partition function is given by 

Z A [J] = N A exp-^Sj(hjj) exp^J p J(-p)D A (p)J(p) (1.11) 

where the source J(p) has the same support as K A (p). N A is the normalization constant 
for which Z A [0] = 1. 

Actually in the case of an infinite lattice N A diverges; to avoid this infrared divergence 
we could work on a finite lattice as long as we deal with the generating functionals; soon 
we will work with the n-point connected Green functions, with n > 1, in which N A , 
related to the vacuum diagrams, does not appear. At that point we could take the limit 
of infinite lattice. For notational semplicity we will work always with an infinite lattice. 

Expanding formally Sj(hjj) in powers of jj, and using the identity 

e " * I p J(-P)D A (P)J(P) j-Q % S ^ J p J(-p)D A (p)J(p) = 
i=l &J{~Pi) 

flttbH) +hD A (p i )—^—}\ rP{p)=DA{p)J{p) = e^i^ DA{p) ^f[<P(p t )U {p)=DA{p)J{p) 
i=i 0( P\ Pi) i=i 

holding when p±, ...,p n are within the support of K A (p), we get 

Z A [J\ = iV A eU ^^^W e f/^^W4e-i^(^ (!))=DAWJ(!)) (L12) 

One has 

Z A [J] = e^[J] (L13) 

where WaJJ] is the functional generator of the connected Green functions. Substituting 
(1.13) in (1.12) one gets 

^[HJ/MWW = iV Ae t/ P ^)^W4i)e^^WU (p)=DA(p)J(p) (1.14) 

The RHS of this equation is the functional generator of the connected Green functions, 
apart from the free contribution D A (p) of the two-point function. In the following we will 
study the amputated connected Green functions, which are amputated of the propagators 
D A {pi) attached to the external legs. The functional generator for these Green functions, 
apart from the 2-point tree- level contribution, will be called Va[0], and it satisfies 

e +sVUd = N A e^p J ^ DA(p)j ^e-^ SlW . (1.15) 
V A [4>] has the Volterra expansion 

00 1 „ n 

V A [</>] =lniV A + £- / ( f>(-k 1 )..4(-k n )(2n)%(J2k i )V A , n (k 1 ,...,k n _ 1 ) 

n= l n - Jk!,...,k„ 1 
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V\,n(ki, k n -i) is symmetric under permutations of k\,...,k n , where k n = — J2i=i hi 
and it is periodic under k^ — > + for any i = 1, n — 1 and for any fi. 

As in [17|, [18[ we will study the flow equation perturbatively in Ti. The amputated 
connected Green functions are defined as series in h, V^ n = J2i>o $V An - 

Each graph contributing to V^ n (ki, fc n -i) is an amputated connected graph with 
internal propagators D A , and with the vertices contained in Si. While in (1.14 — 15) 
we restricted the support of 4>{jp) to be equal to the one of K A , it is clear from the 
Feynman graph representation that the external momenta can be arbitrarily chosen. 
V An (ki, fc n _i) is a C°° function on the Brillouin zone. 

Differentiating equation (1.15) with respect to A one gets the flux equation 

dV A [<j>] 1 f dD A (p) 5 2 V A + 5V A 5V A j + n dlnN A 



dA 2J P dA 1 6<p(-p)6<p(p) 6<f>(-p) 6<p(py OA 

Differentiating (1.16) with respect to 4>(ki), 4>(k n ) and setting to zero one gets 
the flux equations for the amputated connected Green functions, for n > 2 and I > 1, 



A,n 



OA 



(ki, ...,k n -l) = ^ / 9D Qff > V A,nl\(h, K-UP, ~P) 



fln ni-1 

+ E E E ( E kpi)vi l ^(k P1 , fep( m -i))Vygi(A;p ni , k Pn ) (1.17) 

P ni+n 2 =n+2ii+«2=« »=1 

where X)p is the sum over distinct permutations of fci, k n , with k n = — Yli^ 1 %■ 

In the case / = the first term in the RHS is absent, and for I = 0, n = 2 the zero 
solution must be chosen. 

Since the sums in the RHS of (1.17) satisfy U > 0, rti > 2, for given I and n they 
contain a finite number of terms and can be solved iteratively. Indeed in the RHS of 
(1.17) one has V A %, with s < I — 1, or with s = I and E < n — 2. Ordering the vertices in 
such a way that V A „ preceeds ^ ^i < ^ or ^ = ^ anc ^ n i < n 2, it follows that in 
the RHS of (1.17) there appear only vertices preceeding V A n . Therefore (1.17) is a linear 
differential equation in V^. 

V A \ can be determined assigning its value at a particular value of A. At A = A 
equation (1.15) reduces to 

V Ao [(f)}--SM+n\nN Ao . (1.18) 

which provides a set of boundary conditions on V A n at A = Ao- It follows that at A = 
one has 

i 4 »2t/(0 

v >) = « w ; J E |r|r(o) = (O ; <i(o, o, o) = 7 c> 

(o) = , /?W = , i (0) = -g. (1.19) 
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The flux equations for n = 2 and n = 4 determine the three vectors (a™), and 
(7O) in terms of the three vectors (cf**) , z = 1,2 and 3 given in (1.5); in perturbation 
theory the viceversa is also true, since in this case (1-17) is a first-order linear differential 
equation. 

Instead of fixing the boundary conditions for the flux equation (1.17) on all the vertices 
at A = A , one can equivalently impose the renormalization conditions (1.19) on the 
relevant terms, while keeping fixed at A = A the boundary conditions on the irrelevant 
terms. 

The renormalization conditions (1.19) can be expressed in terms of the renormalization 
conditions on the proper vertices T n , using the relations among the proper vertices and 
the connected Green functions. A particular case is = (3^ = 7® = for I > 1; 
in this case the renormalization conditions (1.19) are equivalent to the renormalization 
conditions (1.7). 

From a practical point of view it is simpler to make perturbative computations with 
proper vertices. However the flux equations are simpler for the amputated connected 
Green functions than for the proper vertices, so we have chosen to prove renormalizability 
using the former. A proof of the renormalizability of 4 using the flux equations for the 
proper vertices has been given in [[18] . 



The Green functions Um^^ooV^l can be obtained from Umx -^ooVx^ using (1.15) 
and the semi-group property of the heat kernel, 

e+*v„M = ^0 e | f p s^l D o-DA](p)^j e+ iv A [4>] n 20) 

Following the same logic as in (1.11 — 15), we can compute Um\ - > . oo VQ l } l i n t erms of 
Feynman graphs, in which the propagators are given by /otia ^oo[-Do — Da](p), and the 
vertices are given by limA -^ooVAn- 

Since m 7^ 0, the support of [Dq — D\](p) is bounded; from this fact and from a 

y(l) 



bound on vjp n which we will obtain later one could easily show that, if limA ^ooVj 



exists, then UmA -^oaVon exists. Therefore we can fix the renormalization conditions at 
A = Ai instead of the conditions (1.19): 

VllM = a? ; \ E ^f^(O) = P? ; OO, 0> 0) = J? (1.21) 



with a\ = f3{ = and 7} = —g. The vectors aS l \ /3® and 7^ can be computed 
in function of the vectors a[ , (3^ and 7^ ; order by order in perturbation theory, the 
viceversa is also true. Therefore for m 7^ the renormalization conditions (1.21) are 
equivalent to the renormalization conditions (1.19). For technical reasons we will work 
with the former. For m = the renormalization conditions (1.21) can be imposed, but 
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it is not trivial to prove the existence of the functional Va for A = 0, due to the presence 
of infrared divergences. This issue is studied in |18], [24| using the flow equation. 



1.3 Proof of renormalizability. 

In the following we will use the symbol P(xi, ...,x n ) to indicate each time a different 
polynomial in x%, ...,x n , with non-negative coefficients independent of Ao- 

Indicate formally by d%vjp n (ki, the differentiation of vff n (ki, k n -i) z times 
with respect to k ilfil , A; in _ lAtn _ 1 , where k ijfl , are external momenta. 

Let V(ki, k n -i) be a differentiable function of the momenta k\, k n -\. In the 
following we will use the Taylor formula 

n-1 ! o 

V(k ll ... J k n . 1 ) = V(0,...,0) + J2J2 / dtk i ,-—V(q u ... J q n . 1 )\ qi ^ (1.22) 

Consider a function f a (pi, ■■■,p n -i) which is symmetric under permutations of p±, p n , 
where p n = — Yh=i Pi] the subscript a indicates some collection of indices. Introduce the 
norm 

A = MaX {a} SupM A \fa{Pl, ...,p n -l)\ ( 1 - 2 3) 



,...,n, 



where Ma is the domain satisfying the conditions pf < ^Max(A 2 ,r] 2 ), for any i = 1 
p n = — YJi ~ l Pi] V is a fixed quantity, whose role will become clear later. Since \p n \ > -\pn\ 
for p u belonging to the Brillouin zone, it follows that a bound on ||/||a implies a bound 
on the norm with the form (1.23), with the domain pf < lQMax(A 2 , i] 2 ). While the norm 
(1.23) depends on A , the latter norm does not. 

An useful property of the norm (1.23) is ||/|| A / < ||/||a» for A' < A" < A . 

Let us now discuss more precisely the conditions to be imposed on the term S{ rr in the 
action (1.5). We will assume that the coefficients S® n (kx,...,k n -i) of the perturbative 
Volterra expansion of Si rr be C°° functions on the Brillouin zone, satisfying the bounds 

\5—n—z a 

II^J|a<-j^-P(1b-^) , 1 = 



\5—n—z A 

l^ilniiA < — » 1 > ( L24 ) 
A Ai 



for any fixed 77 and for Ai < A < A . These conditions, together with the vanishing of 
Va° 2 = S\o 1 2> §i ye ^ ne previously mentioned conditions on Si rr . 

Renormalization consists in showing that d%V$ n converges and is rotation-covariant 
in the limit A — * 00. In order to prove renormalizability we will first prove the following 
proposition. 
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For Ai < A < A one has 

\\dld;d z k vi s l\\ K < VA^ M + ^(ln A, (1.25) 

for any s, for i>, w — 0, 1, z > and n > 2. To prove this we will proceed by induction 
in the loop index. 

Equation (1.25) holds for s = 0, as can be seen by the following power counting. 
Consider a tree-level amputated connected graph with n external legs, I internal 
propagators and V n bare vertices with n legs; at tree level one has n > 4. 
The internal propagators satisfy the bound 

\\dXd z k D A (k)\\ < aA- 2 - w ~ z (1.26) 

where a is a positive constant. 

This graph is bounded by const, x A b , with 

6= £(4-n)F fi -2J (1.27) 

n>4 

where we used the bounds (1.24), (1.26) and the fact that const, x A 4_n > ^-^P(ln ^). 
V = J2n Vn is the total number of vertices in the tree graph. Using the identities 

n + 2I=Y,nV n ; V-I=l (1.28) 

ra>4 

it follows that b = 4 — n. Differentiating with respect to y, the only non-vanishing 
contributions are those in which there is at least one irrelevant vertex V n , with n > 4; 
using (1.24), it follows that \\d y Vl 0) n \\ A < ^P(ln^) . 

Using the fact that K A is a C°° function on the Brillouin zone, it is easy to see that 
differentiating the tree-level Green functions with respect to an external momentum, or 
with respect to A, lowers b by one. Therefore the induction hypothesis is true for s = 0. 

Convergence and rotation-covariance of d^V^l for A — > oo follow trivially. 

Differentiating (1.17) with respect to y and k we get the following inequalities, for 
Ai < A < A , 

IIMWgJU <diA IIWvw2IIa + 

ni+n 2 =n+2 h+l 2 =l vi+V2=v zi+Z2+z-j=z 

where d\ and d ZltZ2 ^ Z3 are constants, independent of A and A . 

Assume by induction hypothesis that (1.25) holds for s < I — 1; then we want to show 
that (1.25) holds for s = I. 
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We will proceed in steps; in the first three steps we make a few remarks which are 
used in the following steps. Then we prove (1.25) in a sequence of points, each point 
being established by using the previous ones. This sequence consists in increasing the 
number of legs and going, for fixed number of legs, from irrelevant to relevant. The first 
point consists in proving (1.25) for s = I, n = 2, z > 3, that is for the irrelevant part of 
the two-point function. 

i) For n + z > 5 and n > 2 one has 

/■Ao \5—n—z A 

\mv}Zl\\A< J a dA'\\d A ,d;d* k Vll%, + ^^P(\n^) (1.30) 

where we used the boundary conditions (1.24) on the irrelevant operators. Therefore if 
(1.25) holds for s — l,w — 1 and n + z > 5 , with n > 2, then (1.25) holds for s — I, 
w = and n + z > 5, n > 2. 

ii) In the inequality (1.29), the only terms in the RHS in which there appear vertices 
of order / are those with n\ < n — 2, n<i > 4 or viceversa, since there are no tree- level 
two-point vertices. If (1.25) holds for s = I up to n — 2 external legs, then by the induction 
hypothesis it holds for s — I, with n external legs and w — 1. 

iii) We will use the following inequality, holding for A x < A, 

l^<(o,-,o)| < |^4(o,..,o)|+ r dA'\\d A ,d;d z k vV n \\ A , . 

J Ai 

iv) The case of the two-point function can be treated in the following steps: 

a) using (ii) for n — 2, and observing that there is no term V® 2 in the RHS of (1.29), 
it follows that (1.25) holds for s = I, n = 2 and w — 1; 

b) from (iii) for n = z = 2, v = and A = A , and using the induction hypothesis, it 
follows that the renormalization counterterms satisfy the bound jc^l < P(ln^). Using 
this fact it follows that (1.30) holds also for s = I and n = 2,z > 3. From (a) one gets 
(1.25) for s = /, n = 2, z > 3 and w = 0; 

c) making the Taylor expansion (1.22) for d v y d 2 k vH(k) and using (iii) it follows that, 
using the renormalization conditions (1.21), (a),(b), and the fact that in the Brillouin 
zone one has < | we get (1.25) for s = l,n = z = 2 and w — 0. 

d) making the Taylor expansion (1.22) for 9^dkV A \(k), using the hypercubic symmetry 
and (c) it follows that (1.25) is true for s — I, n — 2, z — 1, w — 0. 

e) making the Taylor expansion (1.22) for dyV A \(k), using the renormalization con- 
ditions (1.21), (a,d) and (iii), it follows that (1.25) is true for s = I, n = 2, z = 0, w = 0. 

Therefore we have shown that the induction hypothesis is true for s — l,n — 2. 

v) Use (i) for n = A, z > 1 to get (1.25) for s = I, n = 4, z > l,w = 0. Using (iii) 
for n = 4, z > in the Taylor expansion of dyV A ^ 4 (ki, ^2,^3) we get, using the boundary 
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conditions (1.21) and the first part of (v), the fact that (1.25) holds for s = I, n = 4, z = 
0,w = 0. Therefore the induction hypothesis holds for s — l,n — 4. 

vi) Assume that (1.25) is true for s = /, 4 < n! < n — 2; then using (i) and (ii) it 
is easy to prove (1.25) for n external legs. Since (1.25) holds for s — l,n — 4, then by 
induction it holds for s = I, n > 4. 

Therefore (1.25) is true for s — I. By induction, it is true for any s. 

From (1.25) it follows that V^ n is bounded for A — > oo; we have now to prove that 
vjp n converges in this limit. 

Define D ZjA (k) = lim Ao ^ 00 dl : D A (k) and V^ n (ki, k n -i) = UmA ^ oo dl : V^ n (k 1 , k n -i). 

We will now prove the existence of V} l An (ki, moreover we will show that it 
is a continuous function of the momenta and satisfies 

^z,A,n(kl, k n -l) = V},Al,n(kl, ■■■■> kn-l) 

+ \ I ^^A'(p)^i h2 (*l,...A-l,P,-p) 
ni—1 

+E^a,A'( E kpi)V^L, ni (k P1: fep^-i))^,^^, k Pn )} (1.31) 
1=1 

for ra + 2 < 4 and 

11 — 1 

+$>^ a ,A'( E ^Vi'^^Pi, ^-i))^^, k Pn )} (1-32) 
1=1 

for n + z > 5, Here and in the following of this section we will use the symbol J2 = 
J2pJ2n 1 +n 2 =n+2J2i 1 +i 2 =iJ2z 1 +z 2 +z [i =z, where J2p is the sum over distinct permutations of 
ki, k n , with k n = — J2i 1 kj. 

The proof of convergence is made by induction in the loop order I, so that we start 
by considering the case / = 0; in the limit A — > oo the only bare vertex corresponds to 
the tree-level term g0 4 ; the propagator converges to Da; a tree-level graph on the lattice 
converges to the corresponding tree-level graph of the continuum g0 4 theory and it is 
easy to check that (1.31 — 32) hold. 

Assume that the induction hypothesis (1.31 — 32) is true for s < I. Then we want to 
show that it is true for s — I. The integrated flux equations are 

d Z k VH(h, kn-i) = d Z k VP un (ki, kn-l) + 

+ \ C dA ' { I P ^'DM^V^Uh, fc„-i,p, -p) 

"1 — 1 

+E^9 A ,ZM E kpiKW^likpi, kp^-i^V^iikpn,, k Pn )} (1.33) 
i=i 
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for n + z < 4 and 

+i £° dA'iJ d^D^d^V^iku k n -!,p, -p) 

+±d?d A ,D A ,(f; k Pi )d z k ^l(k P1 , A^-i))^^^, k Pn )} (1.34) 
i=i 

for n + z > 5. 

We start by proving convergence for n = 2,z > 3. The RHS of (1.34) depends only 
on lower perturbative order in the loop expansion, so that by the induction hypothesis 
its vertices converge, and the limit of the integrand exists. Furthermore, the first term in 
the RHS of (1.34) goes to zero in the limit. From (1.25) it follows that, for rj > Max\ki\, 
the integrand appearing in the RHS of (1.34) is bounded by a A -independent integrable 
function of the integration variables; by the dominated convergence theorem of Lebesgue, 
it follows that (1.32) holds in this case. 

For n = 2, z = 2 analogous considerations can be applied to the integral of the RHS 
of (1.33); the first term in the RHS of (1.33) is Taylor-expanded using (1.22), using the 
renormalization conditions (1.21) and the convergence proof for n = 2,z > 3; it follows 
that d%vjp 2 converges. The proof is repeated for z — 1 and z = in a similar way. For 
n = 4 one realizes that the RHS of (1.34) contains V^ n , with I' = I but only with n' — 2 
and so using the previous result the proof holds again. 

For n > 4, as in the proof of (1.25), a systematic induction in n starts working with 
(1.34). 

Having proven that the vertices converge, let us now show that they are continuous 
in the momenta. In fact, the integrands in the RHS of (1.31) and (1-32) are continuous 
functions of the momenta by induction hypothesis; using (1.25) and (1.26) for suitable r) 
we get a momentum-independent integrable bound which is a function of the integration 
variables, so that a standard theorem leads to the conclusion that the RHS of equations 
(1.31) and (1.32) are continuous functions of the momenta. 

It is now an easy task to show that the limit for A — > oo and the limits involved in 
the derivatives commute, that is 

V?A,n( fc i> •••> k n-i) = dklimA^ooV^lih, fc n _i) (1.35) 
We give the proof for z — 1. From 

Vll(h + Ak, k n -i) - Vl l l(h, fc n _0 = [ kl+Ak dpd p vi l l(p, k 2 , k n -i) (1.36) 

and by another application of the Lebesgue theorem we have 

,. Vilih + Ak, fc n _Q - Vll(h, k n ^) 



13 



I rk!+Ak (l] 

= AkJk dp limA °^°° d P V ^(P, •••> k n-i) (1-37) 

where the integrand in the RHS is a continuous function, by the previous results; taking 
the limit Ak — > one gets immediately the assertion. 

From (1.25) it follows that lirriAo^oodyV^ = 0; from the y-independence of the bound 
(1.25) one could easily repeat a proof analogous to (1.35) of the fact that the two limits 
commute; therefore the renormalized Green functions are independent of Si rr . This 
independence is a manifestation of the large arbitrariness in the discretization of the field 
theory on the lattice. 

Rotation-covariance of lim^^ooV^ is proven inductively in the loop index, using the 
fact that in (1.31 — 32) the coefficient functions D z \ are rotation-covariant. 



We have shown therefore that lim^^^vjl^ki^ & n _i), with the renormalization 
conditions (1.21) are renormalized rotation-covariant Green functions, satisfying the 
bound (1.25), order by order in the loop expansion; furthermore they are independent 
of S irr . For m ^ the Green functions /^a ^oo^2(^i, &n-i), with renormalization 
conditions (1.19), can be obtained from lim\ Q ^ 00 V < ^ n {ki ) fc n -i) using (1.20) and there- 
fore exist, are rotation-covariant and are independent of S irr . This completes the proof 
of renormalizability. 

1.4 Generalizations. 

Let us now relax the condition required previously on the tree-level two-point irrele- 
vant terms. It is easy to take into account irrelevant terms with two legs at tree-level by 
including them in the kinetic term. Add for example the term 

A5[0] = 2Ali 0( ^ )(P)20(A;) (L38) 
to the action (1.4 — 5). The propagator replacing (1.8) is then 

D,(P) = ■ (1-39) 

2 

The parameter £ is chosen to satisfy £ > — j^, in such a way that 

\Ddp)\<-^— 2 (1-40) 

pZ _|_ m z 

where (3 is a positive constant. 

In this case the propagator D^(p) has a single maximum in the Brillouin zone, and it 
coincides with D{p) in the limit A — > oo. 
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The cut-off propagator is now defined as D^(p) = Kji{p)D^{p), replacing (1.9). 
Then we can repeat the proof of renormalization along the lines of subsection 1.3, once 
we replace D\(p) with D^{p). It should be noticed that the bound (1.40) has been 
imposed in order to verify the induction hypothesis (1.25) at s — 0, and that £ plays the 
role of y; for instance from 

||<9 ? AulU < T2 

it follows that (1.25) holds for s = 0, with y replaced by £. Therefore in the continuum 
limit the renormalized Green functions do not depend on £. 

A similar analysis holds for other two-point irrelevant terms. One can prove renor- 
malizability provided (1.40) holds and the propagator has the usual continuum limit H, [| 



A last remark on the irrelevant terms is in order. In the formulation of subsection 1.3, 
Si rr is an assigned function satisfying (1.24). In some application the irrelevant terms are 
instead assigned as a function of the relevant terms, for instance one could be interested 
in an action containing 

1 

C 3 4 +-T2C 3 6 + ... 

To prove renormalizability the procedure is the one used in point (iv. b) in subsection 
1.3. One can easily show that the bound 

l4°|<P(ln^) 

holds. The rest of the induction proof then follows as in the previous subsection. 
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2 Renormalization of matter field theories on a lat- 
tice. 



2.1 Renormalization of matter field theories on a lattice. 

The renormalization of a general class of matter field theories with a finite number of 
physical fields and of auxiliary fields on a space-time lattice (not necessarily hypercubic) 
can be done in a similar way to the one followed in the first section. We will not reproduce 
all the steps followed there; we will deal only with the points which are significantly 
different in the two cases; we will restrict our attention to rectangular lattices with 
lattice spacings a\ > a 2 > a 3 > a 4 . We define p^ = -^-sin^^- and Aq = ^ ^ a~ 2 . 

Consider a lattice theory in which there are fields <3>j with canonical dimension A« = 
1, | or 2, for scalars, fermions and auxiliary fields respectively. For semplicity all fields are 
chosen to be real. It is straightforward to generalize the analysis to the case of complex 
fields. 

The action reads, in momentum space, 

m = \ln m-pW(p)®m + Sim (2.i) 

The matrix propagator is a periodic C°° function on the Brillouin zone. 
Introduce the propagator with infrared cut-off, 

D A:ll (p) = D lJ (p)K A (p) (2.2) 

where K A is defined as after (1.9). 

Assume that the propagator satisfies the following bounds: 

IK^£>A,iil|A < aA A ^- 4 - z ~ w (2.3) 

for < A < A , and a a positive constant depending on z and w. 

Si contains the interaction terms, both relevant and irrelevant. The tree- level part 
of Si is at least cubic in the fields; each relevant tree-level term of Si has a parameter 
which is independent from A . 

The irrelevant terms are chosen to be linear in the parameter y. 

As in the first section, introduce the functional generator V A [Q] for the amputated 
connected Green functions (apart from the tree-level 2-point contribution); each graph 
contributing to these Green functions depends on A only through its internal propagators 
D A . Therefore V\ satisfies a flow equation analogous to (1.16). 

Indicate with V^j, where I = i±, ...,i n , the amputated Green function with external 
legs of fields $ il ,...,$ itl . 
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The amputated connected Green functions are defined as series in h, V A>I = J2i>o ^ Vji} 
The flux equation has the form 

= \jjL 9AD A , jk V^ k + £ d A D Ajk viy uJ Vi l t k (2-4) 

where I± U 1% = I; for semplicity the momenta and the sums have not been written 
explicitly. 

Differentiating the flow equation with respect to y and k we get the following inequal- 
ities, for Ai < A < A , 

IIMWSlU < E^ AAl+Aj_1 II^^WKillA 

+ T,d Z u^ 3 J^A A ^- 5 -^ (2-5) 

where dij and d ZljZ2t z 3 ,h,i2,i,j are constants, independent of A and A . 

Impose the following renormalization conditions on the relevant terms: 

^^(0,...,0)=ag (2.6) 

for TZ =1 A is +z<4. 

These renormalization conditions generalize those in (1.21). If there are symmetries in 
the theory, which can be maintained at quantum level, one can impose relations among 
the renormalization constants present in (2.6); for instance in the case of the scalar 
theory discussed in the first section, the hypercubic symmetry and the discrete symmetry 
4> — > —(f) reduce the conditions (2.6) to the simpler form (1.21). 

Assume that Si„ satisfies the following conditions: 

\\d z A\\A< ^ mn^) ,1 = 

ll^£ >B || A <^— P(ln^) , l>0. (2.7) 

These equations generalize (1.24). 

Let us prove that for Ai < A < A one has 

II^W$IU < VA^^^^ln^ (m^D (2.8) 

for any s, for v, w — 0, 1, z > and n > 1. 

We will prove (2.8) by induction in the loop index. 

The assertion is true for s = 0, as can be seen by the following power counting. 
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Consider a tree level amputated connected graph contributing to V^J. Using the 
bounds (2.3) and observing that each factor A Al coming from an internal propagator 
is cancelled by a factor coming from a vertex, it follows that this graph is bound by 
const, x A 6 , with 

b = 4(F- 7) = 4- (2.9) 

ext ext 

J2ext A« is the sum of the canonical dimensions of the fields appearing in the external legs 
of the graph. 

Differentiating vj^} with respect to y, the only non-vanishing contributions are those 
in which there is at least one irrelevant vertex, satisfying (2.7). Therefore ||9^V^j||a < 
A b+v Av v P{\n^). 

It is easy to see that differentiating the tree level Green functions with respect to an 
external momentum, or with respect to A, lowers b by one. Convergence of d z k V^} for 
A — > oo follows trivially. 

Therefore the induction hypothesis is true for s = 0. 

Assume that (2.8) is true for s < I — 1; then we will show that (2.8) holds for s — I. 
We can follow essentially the same steps as in the first section. Let us only point out 
a few differences. 

The first step is the same, provided n + z > 5 is substituted by J2 \ + z > 5. 

In the second step, in the inequality (2.5) the only terms in the RHS with V^ h are 
those with n — 1 external legs. If (2.8) holds for s = I up to n — 1 external legs, then it 
holds for s = I, with n external legs and w — 1. 

As in the first section, we prove (2.8) by induction in a sequence of points, each point 
being established by using the previous ones. This sequence consists in increasing the 
number of legs and going, for fixed number of legs, from irrelevant to relevant. It is easy 
to prove (2.8) for the tadpoles. The next point consists in proving (2.8) for s = l,n = 2 
and X; A, + z > 5. The rest of the proof is quite similar to the one in the scalar case. 
Analogously one can prove convergence. 

From (2.8) it follows that Um^^oodyV^j = 0; therefore the renormalized Green 
functions are independent of Si rr . 

If the relevant part of Si is hypercubic-invariant and the pointwise limit of D^j 
is 0(4)-covariant in the limit A — > oo, it is possible to choose the renormalization 
conditions in such a way that the renormalized Green functions are rotation covariant in 
the continuum limit. The proof of this fact is along the lines of a similar proof in Section 
1. 
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3 Renormalization of fermionic models. 



In this section we will study the renormalizability of Yukawa theories. We will start 
by reviewing briefly the doubling problem and the Wilson fermions. Since the Wilson 
propagator satisfies the bound (2.3), it is straightforward to prove the renormalizability 
of Yukawa models with this kind of fermions, either using BPHZ 0, |3], || or the Polchinski 
approach here studied. 

Staggered fermions [[11[] do not satisfy the bound (2.3), due to the presence of doublers 



at the border of the Brillouin zone. It has not yet been proven that interacting models 
with staggered fermions are renormalizable. We describe a simple kind of staggered 



fermions, which are staggered in one direction only, and have a single doubler | I2"| , |13 
Although the fermionic propagator does not satisfy the bound (2.3), it satisfies a similar 
bound on a reduced Brillouin zone. Using this fact, we study a chiral sigma model on 
a hypercubic lattice, rewriting it equivalently on a reduced Brillouin zone and showing 
that it is renormalizable as a theory defined on a sublattice. Since renormalizability can 
be achieved maintaining the same translation symmetries as in the theory defined on the 
hypercubic lattice, renormalizability on the sublattice implies renormalizability of the 
chiral sigma model with this kind of staggered fermions on the hypercubic lattice. 

The treatment of the flow equation with auxiliary fields, considered in the previous 
section, plays a crucial role in this proof of renormalizability. 

We conjecture that the renormalizability of Yukawa models with Kogut-Susskind 
fermions can be proven in a similar way by studying them first on a hypercubic sub- 
lattice with lattice spacing 2a, on which the fermionic propagator satisfies the bound 

(2.3) in . 

3.1 Naive fermions, the doubling problem and Wilson fermions. 

The naive fermionic action for massless fermions on a space-time hypercubic lattice 

is 

a 3 4 - 

J o = tt 12 12 ^xi^x+i, - V'x-a) ; (3- 1 ) 

/ X fj,=l 

the gamma matrices are hermitian and satisfy {7^,71/} = 25^ u . 
The inverse propagator for the naive fermionic action is 

S-\ P )= i j2l,P, (3.2) 

where p^ = ^sin^p^a); S' 1 ^) has zeroes for sinp^a = 0, that is for p^ = 0,^; it 
describes 16 Dirac fields in the continuum limit. Iq shares this property with any bilinear 
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and translationally invariant fermionic action, whose propagator satisfies the following 
properties p0| : 



i) reflection (9 ) symmetry: S'~ 1 (pi,P4) = 74>S' _1 ^G i, — ^4)74; 

ii) hypercubic space-time lattice symmetry, i.e. invariance under ~ rotations of the 
coordinate axis; 

hi) chiral symmetry: S 1-1 ^) = — 75>S'~ 1 G )75; 

iv) locality, in the sense that S~ 1 (p) is continuous with its hrst derivatives. 

From (i) and (ii) it follows that S^^p) = j^S -1 ^ (Rpp)^^, where is the reflection 
operator on the fi-th coordinate, (R^x) u = (l — 25^ u )x u . From this fact and (iii) it follows 
that S* -1 ^) = — S' _1 (— p), which together with periodicity p^ = p^ + — gives S" -1 ^) = 
for p M = 0, ^; therefore a propagator satisfying these conditions propagates 16 modes. 

This degeneracy is present also in the naive fermionic action for massive fermions. 

Wilson || eliminated this degeneracy introducing a term which breaks the chiral 
symmetry in a hard way (since it has dimension 5), 

3 4 

I\V = ~7T 12 i'xi^x - Tpx+p. - <Az-a) • (3-3) 

^ X fi=l 

The Wilson action is site-reflection positive for r = 1 and it is link-reflection positive 



for < r < 1 p2l ( for a review see [23] ). It describes one light mode and 15 massive 



modes, with masses of the order of ~, which decouple in the continuum limit. 

The propagator for the Wilson fermion is, in the massive case, 

s( ) _ ~i S„ 1,P, + M(p) 

b[P) ~ p 2 + M 2 (p) (3.4) 

where 

M(p) = M+ —p 2 (3.5) 

and r is the Wilson parameter, satisfying < r < 1. M is the mass of the light fermionic 
mode. We recall that p M = ^sin^rr- 

The fermionic propagator becomes the Dirac propagator in the continuum limit a — > 0, 

lim a ^S{p) = — - (3.6) 

1 Hp. Ta«Pa* + M 

The denominator of S(p) satisfies the bound 

1 r~ 2 

< -^7T (3-T) 



p 2 + M 2 (p) ~ p 2 + M 2 
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where r is the Wilson parameter. To prove this inequality, observe that 
f + M 2 (p) = r 2 (p 2 + M 2 ) + (1 - r 2 )[M 2 + p 2 } 

Y^Pl] + raMp 2 > r 2 (f + M 2 ) 



r 2 a 2 



' - ' ''/-J 



for < r < 1 and M > 0. 

It follows that the Wilson propagator satisfies the bound (2.3). Therefore, according 
to the previous section, one can construct renormalizable models using Wilson fermions. 



3.2 Fermions with minimal doubling. 

On a hypercubic space-time lattice, translation-invariance, locality, chiral symmetry 
and CPO (charge-conjugation x parity x reflection) invariance of the action imply the 
existence of doubling of the Dirac modes on the lattice ||. 

In order to have the minimal doubling allowed under these assumptions, either reflec- 
tion symmetry or hypercubic invariance must be dropped. 



In [12, h3| it has been given an example of lattice translation invariant and chirally 



symmetric fermionic action which breaks the hypercubic space-time symmetry to cubic 



symmetry and which has minimal doubling. The action is [13] 

icc~* \ — — 

I f = I + —— i'xji (2^ x - ipx+fj, - 4>x-fi) + Ma 4 Y ^x^x (3.8) 

^ X p£\ X 

where for M^O the chiral symmetry 

Vfe-e^* ! ^-^ 75 (3.9) 

is softly broken, since ip x ip x has dimension 3; Iq is the naive fermionic action (3.1). 
The inverse propagator is 

on 

S-\p) =iY,PM + ^Phi + M. (3.10) 

where p\ = Y^/j^iPt- For A > 1/2 there are only two propagating modes, around 
p = (0,0,0,0) and p = (^,0,0,0) |T2|; in fact the inverse propagator is equal to M 
provided sinp^a = for /i 1, which means p M = 0, - for \i ^ 1; and provided sinpxa + 
A — cosp M a) = which cannot be satisfied if A > 1/2 and p lt =- for some ju ^ 1. 

This is an example of fermions which are staggered in one direction only. 

For | < A < 1 the action (3.8) is also link-reflection positive, so that one can construct 
a transfer matrix over a Hilbert space of physical states [13| . 



The hypercubic symmetry is broken to the cubic symmetry in the directions £2,^3 
and X4, which includes the axis-inversion symmetry ip x — > i^f^tfpn x , with // ^ 1. 
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The action (3.8) has a discrete symmetry reflecting the fermion in its mirror fermion: 

^.-*(-r iA Vi^ ; $ x ^ {-y i/a $ Rlx . (3.11) 

The action is link-reflection invariant, that is, invariant under the antilinear mapping 

9ipx,t = $x,i-tli ; ®4>x,t = l4ipx,i~t (3.12) 

( it is also site-reflection invariant, but not site-reflection positive ). It is invariant under 
CP transformations 

ipx,t -> 74C^f_ £it ; 4>x, t -> -^i-^C' 1 ^ ( 3 - 13 ) 
where C is the charge conjugation matrix. Therefore the action is CPG-invariant. 

The propagator (3.10) does not satisfy the bound (2.3), since it is equal to M -1 in the 
point p = (-, 0, 0, 0). Similar problems arise with the Kogut-Susskind staggered fermions 



This difficulty can be overcome redefining the fermionic variables ip x and ip x , defined 
on the hypercubic lattice S in terms of new fermionic variables defined on the sublattice 
S' characterized by — even. Define, for — even, 



% — % — — 

$2,x = Ipllbi^x -^x+l) ; ^,x = 2^~^+i) 7l7s ( 3 ' 14 ) 



This redefinition is similar to the one made in |19[ , where the Kogut-Susskind fermions 
are rewritten on a sublattice of lattice spacing 2a. 

The fermionic action (3.8) becomes, in these variables, 

J F [tt, *] = 2a 3 ^'^[^ \l» ® 1(^+A - *x-a) + ® 1(^+21 - **- 2 i) 

i A 
-t7s ® <ri(2* a , - V x+2i - ^._ 2i ) + i- 7i ® ^s(2*x - ^x+a - **-a) + aM^.](3.15) 

where 2 X * s the sum over x e anc ^ where = (\E^ i!B ),i = 1,2. cr^ are the Pauli 
matrices. 

The inverse propagator becomes 

s-^p) = r + M = 

7^ ® + t^7i ® 1 sin(2pia) - 075 <g> axp\ + ^71 ® a 3 p5_} + M (3.16) 
^1 ^ a 2 
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which is periodic on the reduced Brillouin zone \pi\ < ^, < ^, /j, ^ 1. In the 
continuum limit this propagator becomes the standard Dirac propagator. 
Let us show that this propagator satisfies the bound (2.3). Observe that 



\\k a s\\ a = \\k a ^-^\\ a <<a\\k a - 1 



T 2 -M 2 ""- 11 al + E 1 



where 



and 



a 2 



B = a\p 2 ± (l ® a 3 pi cospia — 027175 ® 02P1) 

The smallest eigenvalue of B is /3 = — aXp^pi]. 
Therefore the minimum eigenvalue of al + B is 

,2 



a + > M 2 +pl+ p]_{l - aX^l) + ^-[X 2 (pD 2 - . ( 3 - 17 ) 



^1 

A detailed study of the RHS of (3.17) shows that in the range \ < A < 1 the following 
inequality holds: 

a + P>M 2 + l -{2\-l) 2 {p\+pl) . 

Using the norm (1.23) and the inequalities -^p 2 < p\+ p\ < p 2 , for p belonging to 
the reduced Brillouin zone, it follows that 

||AT A : 1 J U < const, x A' 2 (3.18) 
al + B 

and therefore the propagator satisfies the bound (2.3) for z = w = ; it is easy to see 
that (2.3) is satisfied for any z and w. 

The action (3.15) has the link-reflection symmetry 

*s,i-t74 ® 1 ; %,t -> 74 ® Wx,i-t (3.19) 
It is invariant under the CP transformation 

*x,t -> 74C ® <7 3 \iLx lt ; -> ^ £)t C- 1 74 ® cr 3 (3.20) 
and it has the chiral symmetry (for M — ) 

^ -> e^ 760 * 3 ^ ; ^ -> ^e^ 75 ®' 73 (3.21) 
The discrete symmetry (3.11) takes the form 

V x -> 7i75 ® <T2*fl!x ; -> -^i? ia; 7i75 ® <t 2 (3.22) 
which mixes the ^-inversion symmetry and the flavour symmetry. 
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3.3 Renormalizability of a chiral sigma model. 

Let us describe a sigma model on a hypercubic lattice, with two charged fermions and 
two neutral scalars, with a U(l)v vector symmetry and a U(1)a axial symmetry. The 
tree-level action is 

J (0) = I f [rl>,$] + I B [a] + I B [ir] + ^^(u, + 175^)^ + 9i{o 2 x + n 2 x ) 2 } (3.23) 

X 

where If[ip,ip] is the free fermionic action (3.8) with M = 0; 

2 2 

Of 1 TfL 

J s[0] = — - (px+v - <i>z-ix) + a 4 J2 -ir^l 

is the free bosonic action for <p = <j,tt . a x is a scalar, ir x is a pseudoscalar. 

The action (3.23) has the U(l)y phase symmetry ip x — * e ta ip x , ipx e~ ta ?Px', the 
scalars are invariant under this symmetry. It has an axial symmetry U(1)a', the chiral 
transformation on the fermions is given in (3.9); the scalars transform in the following 

way: 

a x + iir x -> e- 2lf3 (a x + m x ) (3.24) 

The action (3.23) has the discrete symmetry (3.11) on the fermions, and a x — > 
or^x-, t^x ^Rxx on the scalars; it is invariant under CP and under cubic space-time 
rotations and inversions in the X2X?,x± directions. It has the link-reflection symmetry 
defined by (3.12) on the fermions, and by Q<J X}t — > Cx,i-t, 07T T 

,t ^xA-t on the scalars. 

The fermionic propagator does not satisfy the bound (2.3), so that it is not straigh- 
forward to prove the renormalizability of this model. 

As in the previous subsection, we will rewrite this model on an anisotropic sublattice, 
in such a way that the fermion and its doubler become two fermions which are coupled 
by an irrelevant term, with propagator satisfying the bound (2.3). For consistency, it is 
also necessary to rewrite the scalars on this sublattice; a scalar splits in a propagating 
mode and in an auxiliary field. Define, for X\ even, 

1 1 

K = — \o x + a x+i ) ; B x = -{ir x + ir x+i ) 

1 1 

F x = -a'\a x -a x+i ) ; G x = -a -1 ^ - ir x+i ) (3.25) 

while the fermions are defined as in (3.14). 

The free A — F part of the bosonic action Ib[&] becomes, on the sublattice S', 

WtA F] = 2a^{^M^ x - A x+2i - A x _ 2i ) + \f x {QF x + F x+2i + F x _ 2i ) 
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+ 7^-^r(Ar+2i ~~ Ar-2i) + — ^ ^2(A X (2A X — A x+fl — A x _p) 



2 

WF X (2F X - F x+ fi - F X _ A )) + — + a 2 F 2 )} (3.26) 
The interaction term is, at tree level, 

42 = 2a A Y^ {91^ x[(A x + ryi7 5 <g> <7iaF x ) + ry 5 <g> + 27175 <g> aiaG x )]^ x 

X 

+g 2 [(Al + BI + a 2 F 2 + a 2 G 2 x f + 4a 2 (A x F x + ^G,) 2 ]} (3.27) 
The tree-level action (3.23) is equivalent to 



/(°) = I[A, F] + I[B, G] + J F [tt, t] + I 



(o) 



*x - [^ /~ t7 ?" U1 + ^ +2 i ^ t7 ?" U1 ]l ® (3-28) 



The translation transformation by one site in the direction xi on S become, on S', 

4r - ^(Ac+21 + 4r) + |(^+2i - ^) 

r 1 - «7i75 ® <Ti iT , , 1 + 27175 ® o"i . 

2 

1 - 2 7175 ® o-j = 1 + i7i75 ® o"i - 
2 + x+2i 2 

and similarly for B x and G^. All these symmetries are preserved in presence of the cut-off 
K\(p) = K ^ 1 -^ 1 ), introduced for all the propagators. Therefore they are compatible 
with the flow equation. Choosing appropriately the renormalization conditions on the 
relevant terms and the boundary conditions on the irrelevant terms, it follows that the 
connected Green functions have the translation symmetry (3.28) for any A. 

Modulo irrelevant terms, the most general quadratic scalar term in the sector A — F 
of the bare action, which is invariant under the transformation (3.28), has the form 

Iff = 2a 4 E'{ f (Al + a 2 F 2 ) + ^±^A X (2A X - A x+2i - A x _ 2i ) + 2a 2 F 2 

"I ^ F X {2F X — F x+2 i — F x _ 2 i) + —F x (A x+2 i — A x _ 2 i) 

+ 7T2- Y,(M2A X - A x+[1 - A x ^) + a 2 F x {2F x - F x+fl - F x ^))} (3.29) 
Za &i 

Due to the axial symmetry the corresponding term in the B — G sector has the same 
coefficients cci, a 2 , as, a±. 

The same arguments lead to the quadratic fermionic term of the bare action of the 
form 

IlT = 2a 4 E'*x[A E h» ® ~ + W*** - *x-2i) 

-^/^275 ® <Ti(2* a - t> x+2i - * x _ 2i ) + ifoli ® + (3.30) 
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The bare interaction term has the same form as in (3.27), with the coupling constants 
g 1 and g 2 substituted by two bare parameters Xi an d X2- Therefore starting with the 
bare lagrangian the renormalization procedure will require 10 independent conditions, at 
any loop order. 

A parametrization of the two-point solution of the flux equation in the A — F sector 
having the bare action just discussed as boundary condition is 

Vh = A^(p)+ A A (p)cos Pl a ; V A F = a 2 (A A (p) - A A (p) cos Pl a) 

V A A = -ia A A {p) sm Pl a (3.31) 

where A A (p) and A A (p) have period ^ in the directions P2,P3,Pa and are, respec- 
tively, periodic and antiperiodic in the direction Pl , with period -. They correspond 
to \[V A ( Pl , ..) + V A ( Pl + = ..)] and \[V A ( Pl , ..) - V A ( Pl + = ..)], where V A is the two- 
point amputated connected Green function on the lattice S. 

Equation (3.29) can be obtained from (3.31) with a suitable choice of A + and A_ at 
A = A . 

Equation (3.31) leads to the following renormalization conditions 

Vh>) = cf ; f,vh m (0)=4(af + a f) ; -£y£i<»(0) = 2cf 
ap l ap k 

V A f(0) =0 ; V^ (0 (0) = 4af + aV 1 (/) ; A^ (0 (0) = ,k^l 

dpk 

A-V£f\<d) = 2iaf ] (3.32) 

Analogous renormalization conditions are imposed in the B — G sector with the same 
coefficients a'i \ a'} l \ The irrelevant terms in (3.32) are kept to preserve the symmetry 
(3.28). 

The renormalization conditions consistent with (3.9), (3.11 — 13), (3.28) and with the 
axial reflection symmetry assign the following expansion of Vq,^ in the neighborhood of 
p = 0: 

V& (0 (p) = 2rf° 5> 7 * + 2(/3? Pi + i/#°)7i ® ^3 + 0(p 2 ) (3.33) 
k 

For the three-point functions the renormalization conditions are 

1(1) 

y£ l(,) (o) = x2° ; ^ , >(0) = x5° ; ^>(o) = ^- (3.34) 

As the fermionic propagators, also the bosonic propagators associated to (3.26) sat- 
isfy equation (2.3), with A a = A B = 1 and A^ = A G = 2; therefore the model is 
renormalizable on the lattice S'. 
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3.4 Back to the hypercubic lattice. 



We consider now the relation between the renormalization on the cubic lattice S' and 
the hypercubic lattice S. 

This relation has some non-trivial feature already in the case of a pure scalar theory. 
Let us consider a scalar theory on S, which is anisotropic in the X\ direction. Then one 
must impose three independent renormalization conditions on the two-point function at 
every loop order. Considering instead the corresponding theory on the lattice S', the 
scalar field is replaced by a scalar field and an auxiliary field, as in subsection 3.3; one 
must impose four independent renormalization conditions, as in (3.32). On S' there is one 
renormalization condition which is not necessary on S; in fact V^ (0 (0) =Aaf + a 2 af ] 
corresponds to 

<y-,o,o,o) = -| a f + 2a ;« 

This extra renormalization condition, which is a requirement on the behavior of the two- 
point function at the border of the Brillouin zone, is allowed by the fact that the a^-term 
and the a^-term in the action (3.29) are proportional, modulo an irrelevant term which 
on S has the form 

/ 4>{-p)<t>{p)(p\-pl) • (3-35) 

Observe that this functional, which is irrelevant on S, is transformed into a relevant 
operator on S'; in fact the limit pi — > ^ on S corresponds to the limit p\ — > on Ef 
and in this limit this operator does not vanish. Therefore renormalizability of the scalar 
theory on S' is a stronger requirement than renormalizability on S. 

Let us show that for the chiral sigma model the renormalizability on S' implies the 
renormalizability on S. Consider on S' a component of the Volterra expansion of V, 

L- L 6 ^^-*-» v " ■■■ < M6 > 

where $/ stands for the fields A, B, F, G, ^. f B / 2 stands for the integral on Brillouin 
zone corresponding to S'; 5^/2 is the corresponding Dirac delta-function of momentum 
conservation. Using equations (3.14) and (3.25) in momentum space, the amputated 
Green functions on S will be suitable linear combinations of the amputated Green func- 
tions on S'. To show this, consider for example the following manipulation in (3.36), 

~ a dp 1 ...V 1 (p 1 , ...) = 2-3 ~ a ..ty(p u ...)(1 + e- pia ) + ^(Pi - -, ...)(i - e- ipia )] 
Jo Jo a 

= l^dp 1 ...^(p 1 ,...)(l + e-^ a ) (3.37) 
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where we have used the periodicity of ^ for the terms not explicitly written. The integrals 
over the momenta p± is understood. Analogously 




a 



d Pl ...A( Pl , ...) = 2-3 Jl d Pl ...a(Pu + e- tpia ) 



(3.39) 



a 



IJ d Pl ...F( Pl , ...) = 2-3 d Pl ...a-V( Pl , ...)(1 - e-^ a ) 



(3.40) 



a 



and similarly for B, G, it. Using (3.37 — 40) for any pi-integration variable, (3.36) can be 
transformed in a multiple integral over the full Brillouin zone. 

The delta-function of momentum conservation on the reduced Brillouin zone can be 
expressed in terms of the delta-function on the full Brillouin zone, 



This last formula shows that in general, due to the presence of the translated delta- 
function, equations (3.36 — 40) do not lead to the Green functions of a theory on S. 
But in our case, due to the imposition of the symmetry (3.28), this actually happens, 
indeed in the sum all the terms containing the translated delta-function must add to zero. 
Therefore we can conclude that every term (3.36) leads to a contribution to the Green 
functions of the hypercubic lattice S convergent in the continuum limit. 

Moreover, (3.38) shows that, concerning the fermionic lines, not only the 'naive 
limit' a — > of the Green functions on S at fixed momenta exists, but also the limit 
lim a -^oV^„{pi + ^, ..) exists and it is not vanishing; its physical interpretation is of course 
the existence of a second fermionic particle. 

3.5 0(4) invariance. 

In order to get an easy proof of the 0(4) invariance of the theory in the continuum 
limit, a further change of the field variables is suitable. Indeed on S' the pointwise limit 
for a — > of the propagator in the scalar sector is not manifestly 0(4)-invariant, as can 
be seen by Fourier-transforming equation (3.26). 

The amputated Green functions of the field a on S have contributions from both the 
Green functions with A and F external legs; thus only in the sum, in the continuum 
limit, the symmetry will be restored. 



<5s/2(Pl,P2,P3,P4) = 5 B ( Pi , P2 , P 3, P a) + $b{Pi + ~,P2,Ps,P4) 
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Let us make the redefinitions 



A(x) = A'(x) ; F(x)=F>(x)-±-(A' x+2i -A> x _ 2i ) (3.41) 

oCb 

Analogously for the fields B and G. In terms of the new fields equation (3.26) in mo- 
mentum space becomes, up to irrelevant terms, 

I[A', F'] = \ j ^[A\p)A'{-p){p\ +pl + m 2 ) + AF\p)F\-p)] (3.42) 

Moreover F' and G' (as F and G before) appear in the interaction part of the action 
only in the irrelevant terms. From this fact one can easily show that in the continuum 
limit the fields F' and G' decouple, in the sense that, up to the F' — F' and the G'—G' two- 
point terms, all the interaction amputated Green functions with some F' and G' external 
legs vanish. To show this, let us multiply by a parameter y all the interaction terms in 
the bare lagrangian which depend on F' and G'. For y = 1 the action is equivalent to the 
one on the hypercubic lattice S, while for y = this is not true but the fields F' and G' 
trivially decouple. In general, for y ^ 1, renormalization requires the imposition of the 
renormalization conditions on all the terms with dimension lower or equal to 4. However 
we choose to impose again the conditions corresponding to (3.32), (3.33) and (3.34). The 
theory defined in this way is the original one for y — 1; after taking the limit a — > the 
solution of the flow equation does not depend on y, as it was shown in subsection 1.3, 
and in the remark in subsection 1.4; they are in particular equal to those of the y = 
case. 

We showed in subsection 1.3 that in the limit a — > the solutions of the flux equation 
satisfy the system of integral equations (1.31 — 32). Now we have a system of equations 
involving functions with only A', B', \l/ external legs, while the two-point functions of 
the auxiliary fields are momentum- independent. 

Only 0(4)-invariant propagators appear and then, as in subsection 1.3, we conclude 
that with a suitable choice of the constants of the renormalization conditions (for instance 
all vanishing for 1^0) the continuum limit of the solutions of the flux equations on 
H' are 0(4)-covariant. Thus from equations (3.37 — 39), the terms on 5 coming from 
the functions on 5' with A', B', ^ legs will be 0(4)-covariant; from (3.40) Vf>f' could 
contribute only to V acr , but an explicit computation shows that in the continuum limit 
this contribution vanishes. Therefore 0(4)-invariance on S is proven for an arbitrary 
choice of the extra renormalization parameter discussed in the previous subsection. The 
renormalized Green functions do not depend on this parameter. 

In this section we have studied a chiral sigma model, which has massless particles; 
due to the infrared divergences it is by no means trivial to show that the existence of V\ 1 
implies that of the functional V at A = 0. In the formalism of the flux equation (on the 



continuum) this problem has been studied in |TB, 23] . 
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Conclusions. 

In this paper we have shown that the Polchinski approach to renormalization can be 
applied on the lattice. Using this method we have proved the renormalizability of the 
scalar and Yukawa models on the lattice with Wilson fermions. A proof of these facts 
has been obtained previously by Reisz 0, |3|, f| using BPHZ. We have also treated a 
Yukawa model with a simple kind of staggered fermions. We proved that it is renormal- 
izable provided a renormalization condition on an irrelevant term is added to the usual 
renormalization conditions; this extra condition is related to the two-point scalar Green 
function at the point on the border of the Brillouin zone in which there is the doubler 
pole of the fermionic propagator. While we have shown that such a renormalization pro- 
cedure is sufficient to prove the renormalizability, we have not shown that it is necessary. 
Since the corresponding counterterm (3.35) does not appear in the one-loop bare action, 
only from a three-loop computation on the hypercubic lattice one might see if this extra 
renormalization condition is necessary. 

It would be interesting to investigate the renormalizability of models with Kogut- 
Susskind staggered fermions along the same lines on this paper. For instance in the 
case of four flavours (that is in the case of one naive fermion) we expect that a proof of 
renormalizability can be given, involving fifteen extra renormalization conditions, related 
to the two-point scalar Green function at the points on the border of the Brillouin zone 
in which there are the fermionic doublers. 
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